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ABSTRACT 



Finite amplitude acoustic standing waves in a rec- 
tangular air-filled cavity with various wedge-shape 
boundary perturbations were studied both experimentally 
and theoretically. 

The experimental results showed that geometrical per- 
turbations alter the finite amplitude behavior of the 
cavity and that the nature of these changes are in qualita- 
tive agreement with the predictions of the theory. 

However, quantitative agreement was not observed possibly 
because the perturbation chosen did not satisfy all the 
assumptions of theory. 

There does not seem to be any fundamental obstacle 
preventing the choice of a perturbation that would allow 
the theory to be critically tested. 
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I . 



INTRODUCTION 



Coppens and Sanders [1,2] developed a non-linear 
acoustic model with the dissipative term describing the 
viscous and thermal energy losses encountered at the 
walls of a rectangular rigid cavity. 

Several researchers [3,4] have examined the problems 
and experimental results are in excellent agreement with 
the theory except when a degeneracy exists in the cavity. 

The purpose of this research was to examine the 
finite amplitude behavior at the cavity in the region 
for which perturbation effects can be accurately measured 
and to compare the results with the predictions of the 
model . 
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II. BACKGROUND AND THEORY 



In 1975 [2] Coppens and Sanders formulated a pertur- 
bation expansion for the non-linear acoustic wave equation 
with a dissipative term describing the thermal energy loss 
and the non-ideal nature of cavity. 

This model predicts that when a near-degeneracy exists, 
geometrical perturbations provide a mechanism whereby a 
nearly degenerate mode can affect the finite amplitude 
behavior of the cavity. (Most of the following theoretical 
developments were quoted from [2].) 

The non-linear wave equation of the viscous fluid [2] 



is, 






C ' ‘ 2 'p C 

o o o 



(c 2 n 2 + -2-f ) - - 

1 o I— j at* 1 - ‘ p c 2 at 

*0 o 

where X is an operator describing the physical processes 
for absorption and dispersion. The right side term can 
be interpreted as a virtual source distributed within the 
cavity created by the finite-amplitude interaction of the 
linear waves. 

The standing pressure wave in the rectangular cavity 
has the form 



p ~ cos k X cos ky Y cos k Z cos wt 

X Z 

If the cavity is being driven near resonance, the 
non-resonant terms contribute almost nothing with respect 
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to the resonant terms. So the standing wave pressure 
has the form of 



where = MR cos nk X cos nk Y cox n k Z sin (nwt + <J> ) 

P o<V n x y z v n' 

By substituting (2) into (1) , a set of coupled, non-linear, 
transcendental are obtained; 



R (<>„-©)= NM6Q cos 9 [hi R. R . (<j> .-4) .) 

n sin n n n n j = ]_ 3 n-j sin D Y n-j 



n-1 



-cos 



for all n>l, N 




- £ 
j=l 




R. 

3 



rCOS 

*sin ; 



( 4 > 



n+j 



4>j)] 



(3) 



where D indicates the spatial dimension (D = 1, 2, 3). 

If a geometrically perfect rectangular cavity is driven at 
frequencies near the resonance frequency of the (0,1,0) mode, 
A qq = 0 and only the "family members" (o,n,o) will contribute 
to the finite amplitude behavior. 

In this research we were interested in the behavior when 
a degeneracy exists between the 020 and 100 modes with small 
perturbations on the boundary. By using the boundary condi- 
tion, p*n it can be shown that 

P = p 2w + ep' + e 2 p" + . . . 35 p 2w + ep ' (4) 

where p 2w is the pressure of (0,2,0) mode, e is the pertur- 
bation constant defined as the ratio of a geometric 
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perturbation magnitude (A) to the width (L x ) and p is 
(1,0,0) mode component. 

Equation (4) can be expanded 



P = P 



o2o 



-rpr A oo s ioo Pioo e i<2wt + e + T ° 0> (5) 



where S 10Q = Q 1Q0 sin t 10Q x 



•100 



(1 + 



- — , A is the 

oo 



first coefficient of the Fourier transform of the perturbation. 
Equation (5) describes how much energy of the (1,0,0) can be 
transferred to the (0,2,0) mode. 
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III. DEFINITION OF SOME PARAMETERS 



Several parameters were defined for the purpose of 
simplifying the mathematical formulation without loss 
of physical contents. 



A. STRENGTH PARAMETERS (STRM) 

STRPM = MPQ]_ characterizes the strength of the 
finite-amplitude interaction. Since the microphone sensi- 
tivity (Sm obtained with B&K 4220 pistonphone) is known, 

STRPM = M6Q 1 = 7.07 x 10 _3 

c c 

where M - v '' 3 .— - , S = h(l + = 1.2 with = 1.4 

PC^S.„ v v 

M 

P = 1.293 Kg/m^, c = 345 m/sec 

= RMS output voltage of the first harmonic component 



3. FREQUENCY PARAMETER (FP) 



FP = 



2Q (f-nf. ) 
n l 



nf. 



where f is a resonant frequency of (0,1,0) mode, it is the 

Aw 

driving frequency it is assumed that — — <<1. 

n 

FP indicates where the driving frequency is with respect 
to the fundamental resonance frequency of the cavity. As 
| FP | increases, the non-linear effects become smaller. 
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C. HARMONICITY COEFFICIENT (E(n) 



B(n) = '"-f 1 . 

nf 1 


FPn 

2 Qn 



E(n) characterizes how well the modes of a given 
family are tuned. If E(n) <4 | , the corresponding harmonic 
will be strongly excited. 
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IV. THE PRESSURE DISTRIBUTION WITH A WEDGED PERTURBATION 



y 




The pressure distribution of (0,2,0) mode has the form 



P Q20 = P cos y cos ( 2wt + 0 2 ) 
The equation for the perturbed boundary is 



( 1 ) 



3L 



x = L..u- A [ r~ (y-Kn [0 <y--i*> - u ] } W 



'x- ‘ L y - 4 y 



= l [1+e f (y ,z) ] 



x 



3L 



where £ 



= ± , £ (y,z) =*4-(y-|L ) [u (y j^’-uly-Ly)] 



X 



L '- 1 4 y 

y 
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Differentiating Pq 2 q an< ^ f (y,z) with respect to y, gives 
3 P f 



and 



020 '^2itP 2 tt ^ n s 

= ~l — sin l — Y COS ( 2wt+0 2 ) ' 



[0(y -^, -u,y- V J 

■ ^ (y - I V [5ly ' T 1 ’ 



- <s (y-L y ) ] 



(3) 



Equation (3) is substituted into the following Equation 
(4) , 



9P 1 
3x 



x=L 



r r 3f (Y/Z) ^020 , 
x L 3y 3y J x=L 



x 



x 



giving 



3£l 

3x 



2tt 



3L 



= A sin - — y cos (2wt + 6 0 ){U(y — j*~) 



x=L 



- u (y-L y ) 



8ttPL 



where A = 



, 3L 

+ (y— 4 L y) [ ^ (y -— 4 ^) ~ *5 (y -L y) ] ) 

"X 



(4) 



Equation (4) can be expanded as a Fourier series, 



Hr U=l = A cos <2wt + V A m cos r: y 



x 



m=0 



+ I 
m=l 



b 



m 



sin 




(5) 
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Only A q needed to determine the first order correction term 
due to the degerate between (0,2,0) and (1,0,0) modes. 



1 L v 2 tt 3l x 

A q = — / Y sin r y (u (y— ) - U (y-L ) 

y ° y 1 

+ (y - I L y ) [<5(y-§L y ) - 5 (y-L y )]} dy 
2tt 

So the perturbation correction, P' is 



P* = 



2a A. 



L (£-) 
x C o 



Q 



100 



sin t 



100 



cos 



TTX 



X 



i ( 2wt + 



9 2 + T 100 } 



( 6 ) 



The pressure at microphone position, x=L , can be written as 

following with ^ = (|t) , 

C o L y 

P mic = P 020 + £ P' “ Pte t ei(2Wt + 92) + 7 eA o P 100 sin T i00 x 



e i(2wt + 2 + T 100 ) I 



(7) 



The total pressure amplitude at a microphone position will be 



P • I = P/(l+Bcos T. nn ) 2 + (B sinT ir . n ) 2 



mic 



100 J 



■ 100 ' 



( 8 ) 



where B s 



tt £ A o Q 100 Sln T 100 



15 



[In reference [4], the author made a error in formulating 
Equation (2) ; he used y instead of (y-aL ) and got an 
incorrect value of A 0 .] 
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V . APPARATUS 



A. THE RECTANGULAR CAVITY 

The cavity is constructed from 0.75 inch aluminum 
(Fig. 1) . The interior of dimension of the cavity is 
12.002 inches long, 2..5002 inches high with a width that 
can be varied between 5.50 inches to 7 inches in 0.25 inch 
increments. It was sealed by a thin silicon grease layer 
prior to assembly. There are three different microphone 
positions for the purpose of discriminating the different 
modes . 

Figure 2 shows the pressure distribution for several 
of the lower modes. According to Table I (calculated for 
cavity dimensions, 12"x6"x2 . 5" ) , the (0,2,0) and (1,0,0) 
modes are degenerate. Experimentally (0,2,0) mode was 
about 3 Hz higher than (1,0,0) mode. 

The piston is set flush with the bottom of the cavity 
as near to the wall Y=0 as practically possible, and half- 
way between the wall at X=0 and X=L . In this position 
it can excite the (0,n,0) family of modes but it cannot 
excite the (1,0,0) mode. To drive the (1,0,0) mode it is 
necessary to insert an auxilliary driver (ID-30) at 
Position A. A microphone at Position B will sense the 
pressure of all modes. It determines the resonance 
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frequencies and Q's of the degenerate modes. The microphone 
is placed at Position A to sense the (1,0,0) mode without 
contamination from the (0,2,0) and at Position C to sense 
the (0,2,0) mode without contamination from the (1,0,0) 
mode. 

B. BLOCK DIAGRAM (Fig. 4) 

GR 1161-A coherent decade frequency synthesizer produces 
a driving signal which is precise to +0.001 Hz. This 
signal enters a 2120 MB power amplifier which drives the 
shaker. The piston motion was continuously monitored by an 
Endevco Model 2215 Accelerometer mounted within the piston, 
measured on a HP 400D Voltmeter and the output of the 
accelerometer was observed qualitatively on a Model 13D BR HP 
oscilloscope. A Schlumberger spectrum analyzer was used to 
check for harmonic distortion in the piston motion. 

The output of a microphone with matching preamplifier 
B&K 2801 fed into the three devices; (1) HP 400D VTVM for 
measuring overall voltage leve, (2) the spectrum analyzer 
for analyzing the wave component, and (3) two HP 302A wave 
analyzers for measuring the pressure of the first two 
harmonics . 
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Fig. 1. Adjustable Rectangular Cavity 



TABLE I 



Mode 


Frequency (Hz) 


Mode 


Frequency (Hz) 


010 


565 


130 


2039.4 


020 


1131 


140 


2529.0 


030 


1697 


150 


3032.0 


100 


1131 


111 


2980.5 


110 


1205 


101 


2757.5 


120 


1599.5 


Oil 


2926.4 



x 




(b) 



(c) 



(d) 



Fig. 2. 



(a) Cavity Orientation, (b) Pressure erodes for (0,1,0), 
(c) For (0,2,0), (d) For (1,0,0). 
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Fig. 3. Perturbation Configuration 
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Fig. 4. Block Diagram 



VI. DATA COLLECTION PROCEDURE 



Since the resonance frequencies of the cavity was 
observed to drift with time, the system was allowed warm-up 
for at least one hour prior to data collection. The 
harmonic content of the accelerometer output was analyzed 
and the piston was adjusted until the second harmonic of 
the accelerometer output was at least 50 dB below that of 
the fundamental frequency. Figure 5 shows typical results 
for piston ^2^1 t * le region of the finite amplitude 
measurement. For a driving voltage typical of those used 
in the experiment, measurements were limited to those 
frequencies for which V 2 /V^< 1 %. 

There are three measurement procedures for collecting 
data; (1) a pre-run infinitesimal amplitude measurement, 

(2) the finite amplitude measurement, and (3) a post-run 
infinitesimal measurement. 

A. PRE-RUN AND POST-RUN INFINITESIMAL MEASUREMENTS 

During these measurements, the piston was driven so 
that accelerometer output was less than 0.1 v. The resonant 
frequency of the fundamental and the Q factor first few 
(0,n,0) in the modes were determined from 



23 




Figure 5 
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sit i(«i) 




The resonance frequency and Q of (1,0,0) were determined 

/* 

by driving the cavity at port "A:. 

The harmonicity coefficients were then calculated 

from 



f -nf , 

E = , for E 

n nf, ono 



f r- f 2 



, for E 



100 



The times at which f and f T were taken was recorded. 

U. 1j 



B. THE FINITE AMPLITUDE MEASUREMENT 

Throughout this measurement V-^ had to be kept constant 
which required adjustment of the driving voltage applied 
to the piston. The HP 302A wave analyzer was set to the 
driving frequency and the driving voltage adjusted to 
keep constant as the frequency was changed with the 
microphone in Postion B. The frequency was increased 
0.5 Hz steps through +6 Hz about the resonance frequency 
of a fundamental mode. At each driving frequency, the 
harmonic content of the microphone onput, usually up to 
the fourth, was measured with HP 302A wave analyzer set 
on AFC mode and the time of that measurement recording 
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presented. The results were presented as V^/V^ vs. the 
frequency parameter for a given strength parameter and 
perturbation. 

To determine the exact value of the frequency parameter, 

the value of f at the instant of the measurement of V was 
r n 

determined by interpolation of the values of f found in 
the pre- and post-run procedures. Figure 6 shows the 
drift in f is approximately linear in time. 
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VII. RESULTS 



Figure 7 shows the excellent agreement between theory 
and experiment obtained when there are no perturbations so 
that only one family of modes is excited. 

Some discrepancies were always observed when the driving 
frequency was far from the fundamental resonance frequency 
(for ] FP I < 2 . ) In this region, the piston had to be driven 
hard to keep constant. This causes a second harmonic to 
appear in the piston waveform (as shown by Fig. 3) , thereby 
introducing a linearly-generated second harmonic into the 
cavity. 

For the boundary perturbation of Fig. 3, the results 
are plotted in Figs. 8 through 15. The thin line is the 
theory prediction for the perturbation and the thick line 
includes the perturbation correction. The small circles 
are experimental results. Along the edge of each figure 
is the information about the cavity configuration. 

A. THE WEDGE AT THE CORNER OF THE CAVITY 

Figures 8 through 11 show the results for a wedge at 
the corner of the cavity. The wedge used are those of 
Fig. 3 (a) -(d) respectively. 

B. THE WEDGE AT THE CENTER OF THE LONG WALL 

Figures 12 through 15 show the results for the same 
wedges (Fig. 3 (a) -(h)) but now located at the center of 
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0.177 with 



the long wall. For Figures 12 and 13, A q = 

DELTA for Figure 13 twice that for Figure 12. 

and 15, A =0.500 with DELTA the same as for 
o 

13. 



For Figures 14 
Figures 12 and 
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VIII. CONCLUSIONS 



A. The experimental apparatus and procedures are 
capable of providing data sufficiently precise to verify 
the prediction of the theory in the absence of geometrical 
perturbation if the range of driving frequencies is 
restricted so that the absolute value of the frequency 
parameter is less than 2. 

B. To verify the correction for a geometrical per- 
turbation, it is useful to have: (1) a sufficiently large 

correction to the unperturbed prediction and (2) the sign 
of the correct must change for an absolute value of the 
frequency parameter less than 2. 

C. The experimental results show that geometrical 
perturbation alters the finite amplitude behavior of the 
cavity, and that the nature these change are in qualita- 
tive agreement with the predictions of the theory. However, 
quantitative agreement was not observed. 

D. Sources of the difficulties in obtaining good 
agreement might be; (1) the inability to experimentally 
satisfy conditions of B above, (2) the DELTA function 
arising from the perturbations used may cause trouble in 
the theory. (Also, they may be difficult to accurately 
reproduce experimentally.), and (3) the higher order 
terms neglected in developing the theory may not all be 
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small. For example, terms of order higher than first 

order in Equation (4) of Section II and terms of order 

higher than A_ and all b in Equation (5) of Section IV. 
o m 

E. The theory for the perturbed cavity can probably 
be critically tested, if a perturbation shape is used 
that avoids the difficulties of D(2) and satisfies the 
condition of B. At this time it is not obvious that 
such a perturbation can be found, but there seems to be 
no fundamental obstacle to find a suitable perturbation. 
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Fig. 11 
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Fig. 12 
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Fig. 13 
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Fig. 14 
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Fig. 15 
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TABLE I (Fig. 7) 



mode 


time 




f 

u 


f L 




f 

r 


Q 


E 




FP 


010 

020 

030 

100 

V l = 


0 

3 

6 

9 

lldB 


565.77 

1131.84 

1698.27 

1131.87 

time f 


562.88 564.33 

1127.79 1129.81 

1693.52 1695.90 

1127.59 1129.73 

V 2 V 3 V 4 


195.40 

278.69 

356.51 

264.20 

P 2/ P l 


0 

0.10x10 

0.17x10 

-0.74x10 

Vh 


-2 

-2 

-4 

P 4/ P . 






12 


558.5 


-50.0 


-70.5 


-70.5 


0.0113 


0.0011 


0.0011 


-4.19 






13 


559.0 


-47.9 


-72.0 


-69.0 


0.0143 


0.0009 


0.0013 


-3.85 






14 


559.5 


-46.9 


-73.8 


-67.8 


0.0160 


— 


0.0014 


-3.50 






15 


560.0 


-46.0 


-72.6 


-67.5 


0.0119 




0.0015 


-3.15 






16 


560.5 


-44.8 


-71.2 


-69.7 


0.0204 


0.0011 


0.0012 


-2.81 






17 


561.0 


-43.7 


-69.7 


-73.5 


0.0233 


0.0012 


— 


-2.47 






18 


561.5 


-42.4 


-67.8 


-76.5 


0.0269 


0.0014 


— 


-2.12 






19 


562.0 


-41.0 


-65.7 


-75.9 


0.0316 


0.0018 


— 


-1.78 






20 


562.5 


-39.4 


-63.2 


-82.5 


0.0380 


0.0025 




-1.43 






21 


565.0 


-37.8 


-60.3 


-85.0 


0.0457 


0.0034 


— 


-1.09 






22 


563.5 


-35.9 


-56.6 


-74.9 


0.0569 


0.0052 




-0.74 






23 


564.0 


-34.0 


-52.2 


-67.5 


0.0708 


0.0007 


0.0015 


-0.40 






24 


564.5 


-32.2 


-47.6 


-60.6 


0.0876 


0.0148 


0.0033 


-0.05 






25 


565.0 


-31.1 


-45.0 


-54.5 


0.0994 


0.0200 


0.0067 


0.29 






26 


565.5 


-31.6 


-43.9 


-52.9 


0.0933 


0.0226 


0.0081 


0.64 






27 


566.0 


-33.3 


-46.7 


-56.0 


0.0772 


0.0164 


0.0056 


0.98 






28 


566.5 


-35.3 


-51.0 


-63.2 


0.0610 


-0.0100 


0.0025 


1.33 






29 


567.0 


-37.6 


-55.5 


-70.8 


0.0484 


0.0060 


0.0010 


1.67 






30 


567.5 


-39.8 


-60.0 


- 73.6 


0.0365 


0.0036 


— 


2.02 






31 


568.0 


-42.1 


-64.3 


-74.7 


0.0279 


0.0022 




2.36 






32 


568.5 


-44.1 


-68.5 


- 74.2 


0.0221 


0.0013 


— 


2.71 






33 


569.0 


-46.2 


-69.9 


— 


0.0174 


0.0011 


— 


3.05 






34 


569.5 


-47.2 


-68.7 


-67.0 


0.0156 


0.0013 


0.0016 


3.40 






35 


570.0 


-48.0 


-65.4 


-67.0 


0.0141 


0.0019 


0.0010 


3.14 






36 


570.5 


-48.9 


-65.1 


-67.2 


0.0127 


0.0020 


0.0016 


4.09 






37 


571.0 


-49.8 


-65.2 


-67.5 


0.0115 


0.0019 


0.0015 


4.43 






38 


571.5 


-50.8 


-65.4 


-68.2 


0.0103 


0.0019 


0.0014 


4.78 






39 


572.0 


-51.7 


-65.8 


-68.6 


0.0093 


0.0018 


0.0013 


5.12 






40 


572.5 


-52.5 


-66.3 


-67.4 


0.0084 


0.0017 


0.0015 


3.47 






41 


573.0 


-54.7 


-65.1 


-66.7 


0.0065 


0.0020 


0.0016 


5.81 



mode 


time 


f 

u 


f L 


f 

r 


Q 


E 


010 


44 


566.28 


563.37 


564.82 


194.03 


0 


020 


47 


1132.74 


1128.67 


1130.71 


277.82 


0.94xl0 -3 


030 


50 


1699.77 


1694.97 


1697.37 


353.84 


0 . 17x1 0 -2 


100 


53 


1132.80 


1128.48 


1130.64 


262.03 


-0. 61xl0“ 4 
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TABLE II (Fig. 8) 



mode 


time 


f 

u 


f L 


f 

r 




Q 


E 




FP 


010 

020 

030 

100 

V^-lldB 


0 

3 

6 

9 

time 


568.27 

1132.37 

1701.89 

1132.44 

f 


565.28 

1127.62 

1696.79 

1127.71 

V 2 


566. 

1129. 

1699. 

1130. 

V 3 


77 189.81 

99 238.25 

34 333.14 

07 238.76 

V 4 V P 1 


0 

-0.31x10' 

-0.57x10- 

0.69x10' 

p 3 /Pi 


-2 

-3 

-4 

V*i 




12 


561.3 


-44.3 


-69.9 


-71.3 


0.0216 


0.0011 


— 


-3.95 




13 


561.8 


-43.2 


-69.2 


-72.2 


0.0247 


0.0012 


— 


-3.61 




14 


562.3 


-41.7 


-68.3 


-74.0 


0.0293 


0.0014 


— 


-3.26 




15 


562.8 


-40.2 


-66.8 


-80.0 


0.0347 


0.0016 


— 


-2.92 




16 


563.3 


-38.5 


-65.2 


— 


0.0422 


0.0020 


— 


-2.58 




17 


563.8 


-36.4 


-62.9 




0.0537 


0.0025 




-2.24 




18 


564.3 


-34.3 


-60.3 


— 


0.0684 


0.0034 




-1.90 




19 


564.8 


-32.8 


-57.2 




0.0813 


0.0049 


— 


-1.56 




20 


565.3 


-32.5 


-55.0 


- 72.0 


0.0841 


0.0063 


0.0009 


-1.22 




21 


565.8 


-33.5 


-53.5 


- 67.8 


0.0750 


0.0075 


0.0014 


-0.87 




22 


566.3 


-35.1 


-51.8 


-62.5 


0.0624 


0.0092 


0.0024 


-0.53 




23 


566.8 


-36.9 


-50.9 


-60.9 


0.0507 


0.1002 


0.0032 


-0.19 




24 


567.3 


-37.9 


-51.2 


-60.8 


0.0452 


0.0098 


0.0032 


0.15 




25 


567.8 


-38.6 


-52.8 


-64.6 


0.0419 


0.0082 


0.0021 


0.49 




26 


568.3 


-39.2 


-55.3 


-69.2 


0.0391 


0.0061 


0.0072 


0.83 




27 


568.8 


-40.3 


-58.7 


-75.0 


0.0343 


0.0041 


— 


1.78 




28 


569.3 


-41.5 


-61.7 


-78.5 


0.0300 


0.0029 


— 


1.52 




29 


569.8 


-42.4 


-64.3 


— 


0.0271 


0.0022 


— 


1.86 




30 


570.3 


-44.3 


-67.3 


— 


0.0218 


0.0015 


— 


2.20 




31 


570.8 


-45.5 


-70.8 


— 


0.0188 


0.0010 


— 


2.54 




32 


571.3 


-46.9 


-75.0 


— 


0.0160 


0.0006 


— 


2.88 




33 


571.8 


-48.1 


-76.9 


— 


0.0140 


0.0005 


— 


3.22 




34 


572.3 


-49.3 


-75.2 


-69.1 


0.0122 


0.0006 , 


0.0012 


3.57 




35 


572.8 


-50.6 


-73.5 


-68.0 


0.0105 


0.0007 


0.0014 


3.91 




36 


573.3 


-51.0 


-71.7 


-66.5 


0.0100 


0.0009 


0.0017 


4.25 




37 


573.8 


-51.8 


-71.1 


-66.0 


0.0092 


0.0010 


0.0018 


4.59 




38 


574.3 


-52.5 


-70.3 


-65.7 


0.0084 


0.0011 


0.0019 


4.93 




39 


574.8 


053.8 


-72.3 


-66.7 


0.0073 


0.009 


0.0016 


5.27 



mode 


time 


f 

u 


f 2 


f 

r 


Q 


E 


010 


42 


568.82 


565.95 


567.38 


197.49 


0 


020 


45 


1133.50 


1128.70 


1131.10 


235.40 


-0. 32xl0 -2 


030 


48 


1703.98 


1698.69 


1701.34 


321.25 


-0.48xl0 -3 


100 


51 


1133.85 


1129.10 


1131.47 


238.00 


0.33xl0~ 3 
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TABLE III (Fig. 9) 



mode 


time 


f 

u 


f 2 


f 

r 


Q 


E 


010 


0 


567.85 


565.01 


566.43 


199.52 


0 


020 


3 


1128.61 


1124.81 


1126.71 


296.58 


-0.54xl0 -2 


030 


6 


1697.95 


1693.16 


1695.39 


354.31 


-0.23xl0 -2 


100 


9 


1128.77 


1124.87 


1126.82 


288.63 


-0.96xl0~ 4 


Vj=-lldB time 


f Y 2 


V 3 


V 4 VL V*! 


V P 1 ] 



12 


561.0 


-37.6 


-66.2 


13 


561.5 


-36.2 


-64.8 


14 


562.0 


-34.3 


-62.8 


15 


562.5 


-32.2 


-60.3 


16 


563.0 


-30.2 


-57.5 


17 


563.5 


-29.5 


-55.6 


18 


564.0 


-30.8 


-55.2 


19 


564.5 


-33.3 


-55.5 


20 


565.0 


-35.5 


-55.8 


21 


565.5 


-37.3 


-58.0 


22 


566.0 


-38.9 


-60.5 


23 


566.5 


-40.3 


-61.5 


24 


567.0 


-41.7 


-61.7 


25 


567.5 


-42.5 


-61.3 


26 


568.0 


-42.6 


-60.8 


27 


568.5 


-43.1 


-62.2 


28 


569.0 


-43.3 


-63.7 


29 


569.5 


-43.8 


-65.9 


30 


570.0 


-44.5 


-68.0 


31 


570.5 


-45.2 


-69.6 


32 


571.0 


-45.6 


-70.9 


33 


571.5 


-46.2 


-71.7 


34 


572.0 


-46.8 


-72.7 


35 


572.5 


-47.5 


-73.3 


36 


573.0 


-47.5 


-73.5 


37 


573.5 


-48.2 


-77.0 



mode 


time 


f 

u 


f L 


010 


40 


567.99 


565.11 


020 


43 


1128.94 


1125.12 


030 


46 


1698.43 


1693.67 


100 


49 


1129.10 


1125.21 



— ' — — — — 


0.0468 

0.0550 


0.0017 

0.0021 




-3.84 

-3.49 






0.0684 


0.0026 


-3.14 






0.0876 


0.0034 


-2.79 






- 76.0 


0.1096 


0.0047 


0.0006 


-2.44 


-73.4 


0.1189 


0.0059 


0.0008 


-2.09 


-74.0 


0.1023 


0.0062 


0.0007 


-1.74 


-74.7 


0.0772 


0.0060 


0.0007 


-1.39 


-71.8 


0.0599 


0.0058 


0.0009 


-1.04 


-70.2 


0.0484 


0.0045 


0.0011 


-0.67 


-71.4 


0.0403 


0.0033 


0.0010 


-0.34 


-76.0 


0.0345 


0.0030 


0.0006 


0.00 


on n 


0.0293 


0.0029 




0.36 


ou • u 
-81.3 




0.0266 


0.0031 


0.71 




-80.3 


0.0263 


0.0032 


1.06 






0.0250 


0.0028 


1.41 






0.0243 


0.0023 


1.76 






0.0230 


0.0018 


2.11 






0.0213 


0.0014 


2.46 






0.0196 


0.0012 


2.81 






0.0186 


0.0010 


3.16 






0.0174 


0.0009 


3.51 






0.0162 




3.86 








0.0150 


4.21 








0.0150 


4.56 








0.0139 


4.91 








f 

r 


Q 


] 


E 




566.55 


196. 


79 


0 




1127.03 


295. 


42 -0 


.54xl0~ 2 




1696.05 


355. 


76 -0 


.21xl0 -2 




1127.15 


289. 


76 0 


.llxlO -3 
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TABLE IV (Fig. 10) 



mode 


time 


f 

u 




f L 


f 

r 


Q 


E 




010 




0 


566.19 


563.20 


564.70 


189.30 


0 




020 




3 


1128.86 


1124.95 


1126.91 


287.70 


-0.2x10“ 


■2 


030 




6 


1696.18 


1691.38 


1693.78 


353.24 


-0.18x10 


f 3 


100 




9 


1128.46 


1124.62 


1126.54 


293.29 


-0.32x10 


f 3 


V 1 =-lldB 


time f 


V 2 


V 3 


V 4 


P 2 /P l 


V*i 


Vi 


EP 




12 


558.7 


-42.8 


-69.8 




0.0257 


0.0011 




-4.04 








13 


559.2 


-41.8 


-68.7 


0.0288 


0.0013 


-3.71 








14 


559.7 


-40.8 


-67.9 


0.0324 


0.0014 


-3.37 








15 


560.2 


-39.7 


-66.9 


0.0367 


0.0016 


-3.04 








16 


560.7 


-38.3 


-65.5 


0.0434 


0.0019 


-2.71 








17 


561.2 


-36.7 


-63.5 


0.0519 


0.0024 


-2.38 








18 


561.7 


-34.9 


-61.0 


0.0638 


0.0032 


-2.04 




-76.3 




19 


562.2 


-33.0 


-57.8 


0.9794 


0.0046 


-1.71 








20 


562.7 


-30.7 


-53.7 


-70.7 


0.1041 


0.0074 


0.0010 


-1.38 




21 


563.2 


-28.9 


-50.2 


-65.3 


0.1274 


0.0110 


0.0019 


-1.04 




22 


563.7 


-28.5 


-47.5 


-61.4 


0.1334 


0.0150 


0.0030 


-0.71 




23 


564.2 


-30.4 


047.2 


-59.0 


0.1078 


0.0155 


0.0040 


-0.38 




24 


564.7 


-32.2 


-48.5 


-58.5 


0.0873 


0.0133 


0.0042 


-0. 05 




25 


565.2 


-33.7 


-52.2 


-62.7 


0.0733 


0.0087 


0.0026 


0.29 




26 


565.7 


-35.6 


-56.4 


-70.07 


' 0.0592 


0.0054 


0.0010 


0.62 






566.2 


-37.0 


_ co A 


-74.7 


0.0501 


0.0043 




0.95 




Z / 
9Q 


"JO. He 
_CTQ n 




566.7 


-39.6 


- 77.2 


0.0314 


0.0037 


1.29 




Zo 

29 


-61.5 




567.2 


-42.1 




0.0279 


0.0030 


1.62 








30 


567.7 


-43.7 


-63.7 


0.0232 


0.0023 


1.95 








31 


568.2 


-44.5 


-65.5 


0.0221 


0.0019 


2.29 








32 


568.7 


-45.1 


-66.9 


0.0197 


0.0016 


2.62 








33 


569.2 


-45.6 


-68.0 


0.0186 


0.0014 


2.95 








34 


569.7 


-46.1 


-69.0 


0.0176 


0.0013 


3.29 








35 


570.2 


-46.6 


-69.8 


0.0167 


0.0012 


3.62 








36 


570.7 


-47.0 


-70.4 


0.0158 


0.0011 


3.95 








37 


571.2 


- 47.2 


-70.6 


0.0155 


0.0011 


4.28 








38 


571.7 


-47.7 


-71.3 


0.0146 


0.0010 


4.62 








39 


572.2 


-48.6 


-63.9 


0.0132 


0.0023 


4.95 








mode time 


f 

u 


f I 


i 


f 

r 


Q 




E 




010 


42 


566.35 


563. 


33 


564.84 


186.79 




0 




020 


45 


1129.43 


1125. 


52 


1127.48 


1287.84 


-0, 


,19xl0 -2 




030 


48 


1697.04 


1692. 


26 


1694.65 


353.94 


-0, 


. 78xl0 -4 




100 


51 


1129.12 


1125. 


28 


1127.20 


294.16 


-0. 


. 24xl0~ 3 
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TABLE V (Fig. 11) 



mode 


time 


f 

u 


f L 


f 

r 


Q 


E 


010 


0 


566.90 


563.86 


565.38 


185.80 


0 

-2 


020 


3 


1127.19 


1123.28 


1125.23 


288.08 


-0.49x10 


030 


6 


1695.31 


1690.47 


1692.89 


349.27 


-0 . 19xl0~ 2 

-3 


100 


9 


1126.76 


1122.90 


1124.83 


291.41 


-0.36x10 


V^-lldB 


time 


f V 2 


V 3 


V 4 


P 2 y/P l P 3 //P l 


V*! 



12 

13 


559.4 

559.9 


-39.5 

-37.9 


-66.9 

-65.3 




0.0376 

0.045 


0.0016 

0.0019 




-3.99 

-3.66 






14 


560.4 


-36.4 


-63.9 


0.054 


0.0023 


-3.33 






15 


560.9 


-34.5 


-61.8 


0.0672 


0.0029 


-3.00 


-78.0 




16 


561.4 


-32.4 


-59.1 


0.0856 


0.0039 


-2.66 




17 


561.9 


-30.3 


-55.7 


-72.7 


0.109 


0.0058 


0.0008 


-2.33 


18 


562.4 


-28.7 


-55.6 


-68.2 


0.131 


0.0059 


0.0014 


-2.00 


19 


562.9 


-28.9 


-50.8 


-65.3 


0.127 


0.0103 


0.0017 


-1.67 


20 


563.4 


-30.9 


-50.0 


-62.5 


0.101 


0.0112 


0.0027 


-1.34 


21 


563.9 


-32.3 


-49.7 


-60.3 


0.081 


0.0116 


0.0034 


-1.01 


22 


564.4 


-34.7 


-50.9 


-61.4 


0.065 


0.0101 


0.0030 


-0.67 


23 


564.9 


-36.3 


-54.3 


-67.0 


0.055 


0.0068 


0.0016 


-0.34 


OA 


cr /r cr a 


-37.5 


-59.3 




0.048 


0.0038 




-0.01 


25 


OuD • 4 

565.9 


/ j • Z 




-38.4 


-65.0 


0.043 


0.0020 


-0.32 






26 


566.4 


-39.3 


-69.7 


0.039 


0.0012 


-0.65 






27 


566.9 


-39.7 


- 72.2 


0.037 


0.0009 


0.98 






28 


567.4 


-40.6 


-74.5 


0.033 




1.32 








29 


567.9 


-40.7 


-75.4 


0.030 


1.65 








30 


568.4 


-41.6 


-73.5 


0.030 


1.98 




0.0009 




31 


568.9 


-44.0 


-71.5 


0.022 


2.31 






32 


569.4 


-47.3 


-70.9 


0.015 


0.0010 


2.64 






33 


569.9 


-49.1 


-71.1 


0.012 


0.0010 


2.97 






34 


570.4 


-49.9 


-71.6 


0.011 




3.31 








35 


570.9 


-50.0 


-71.8 


0.011 


3.64 








36 


571.4 


-50.3 


- 72.3 


0.011 


3.97 








37 


571.9 


-50.4 


-72.4 


0.011 


4.30 








38 


572.4 


-50.4 


-72.4 


0.011 


4.63 








39 


572.9 


-51.1 


-71.6 




0.010 


— 


— 


4.96 



mode 


time 


f 


f T 


f 


Q 


E 


u 


L 


r 






010 


42 


566.95 


563.96 


565.46 


189.24 


0 


020 


45 


1127.31 


1123.40 


1125.37 


288.18 


-0.49xl0~ 2 


030 


48 


1695.54 


1690.73 


1093.14 


351.57 


-0.27xl0 -3 


100 


51 


1126.98 


1123.12 


1125.05 


291.39 


-0.27xl0 -3 
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TABLE VI (Fig. 12) 



mode 


time 


f 

u 




f r 


0 


E 


010 


0 


564.94 


562.04 


563.49 


194.44 


0 


020 


3 


1137.75 


1133.55 


1135.65 


270.78 


0.77xl0 -2 


030 


6 


1699.31 


1694.34 


1696.82 


341.28 


0.28xl0~ 2 


100 


9 


1137.85 


1133.56 


1135.71 


264.96 


0.48xl0 -4 



V^-lldB 


time 


f 


V 2 


V 3 


V 4 




P /P 

y i 


i — 1 
0^ 


FP 




12 


557.5 


-53.7 


- 73.9 


- 76.2 


0.0073 


0.0007 




-4.30 




13 


558.0 


-53.3 


- 74.2 


- 76.3 


0.0077 






-3.96 




14 


558.5 


-52.6 


-74.6 


- 77.3 


0.0083 


— 


— 


-3.62 




15 


559.0 


-52.0 


-75.1 


-80.2 


0.0089 




— 


-3.27 




16 


559.5 


-51.2 


-85.3 


-85.0 


0.0098 






-2.93 




17 


560.0 


-50.4 


- 75.3 


— 


0.0107 


— 


— 


-2.58 




18 


560.5 


-49.6 


-74.9 




0.0117 


— 




-2.24 




19 


561.0 


-49.0 


-73.3 


— 


0.0126 


0.0008 


— 


-1.89 




20 


561.5 


-48.2 


-71.6 


— 


0.0139 


0.0009 


— 


-1.55 




21 


562.0 


-47.7 


-69.7 




0.0146 


0.0012 


— 


-1.21 




22 


562.5 


-47.1 


-67.2 


— 


0.0158 


0.0015 


— 


-0.86 




23 


563.0 


-46.3 


-64.7 


— 


0.0172 


0.0021 


— 


-0.52 




24 


563.5 


-45.5 


-61.3 


— 


0.0188 


0.0031 


— 


-0.17 




25 


564.0 


-44.4 


-58.1 


— 


0.0215 


0.0044 


— 


0.17 




26 
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— 
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27 


565.0 


-40.5 
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— 
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— 
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32 
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0.0813 
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35 
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0.0603 
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3.61 




36 


569.5 


-37.8 


-67.0 


-59.7 


0.0457 


0.0016 


0.0037 


3.96 




37 


570.0 


-40.2 


—66.6 


-60.0 


0.0349 


0.0017 


0.0035 


4.30 




38 


570.5 


-42.0 


-66.2 


-60.3 


0.0282 


0.0017 


0.0034 


4.65 




39 


571.0 


-43.6 


-65.6 


-60.5 


0.0236 


0.0019 


0.0033 


4.99 




40 


571.5 


-46.0 


-68.6 


-62.0 


0.0178 


0.0013 


0.0028 


5.34 
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46 
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0.78xl0 -2 




030 


49 
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0. 37xl0 -2 




100 


52 


1139.00 


1134. 


76 1136.88 


268.20 
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TABLE VII (Fig. 13) 



node 


time 
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u 


f L 


f 

r 


Q 


E 


010 


0 


561.56 


558.62 


560.09 


190.57 


0 


020 


3 


1137.45 


1133.15 


1135.30 
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0.14xl0 _1 


030 


6 
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9 
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1133.26 


1135.42 


263.01 


O.llxlO -3 


V 1 =-llc3B 


time 


f 


V 2 


V 3 V 4 


V P 1 


p 3 /Pl P 4 /P ; 



12 
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-53.6 


-74.7 


0.0074 


0.0007 




-4.21 


13 


554.6 
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-74.7 


0.0079 


0.0007 
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555.1 


-52.7 
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0.0007 
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555.6 
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— 
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— 
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-75.1 


0.0114 


0.0006 
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30 
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31 
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— 


2.26 


32 
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0.0047 


— 


2.61 


33 
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-60.9 


0.0309 


0.0032 
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34 
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3.29 


35 
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0.0014 
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36 


566.1 
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0.0010 


— 


3.97 


37 
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— 


4.31 


38 
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-33.2 


-71.0 
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— 


4.65 


39 


567.6 


-31.5 
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0.0017 


— 
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TABLE VIII (Fig. 14) 
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-50.0 
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559.5 


-48.4 
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-47.7 
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560.5 


-47.0 
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561.0 


-46.3 
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561.5 


-45.5 


21 
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-44.6 
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562.5 


-43.6 


23 


563.0 


-42.5 
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563.5 


-41.5 


25 


564.0 


-40.6 


26 


564.5 


-40.0 


27 


565.0 


-40.8 


28 


565.5 


-42.0 


29 


566.0 


-42.6 


30 


566.5 


-41.2 


31 


567.0 


-39.3 


32 


567.5 


-37.0 


33 


568.0 


-35.0 


34 


568.5 


-33.0 


35 


569.0 


-32.6 


36 


569.5 


-33.6 


37 


570.0 


-35.5 


38 


570.5 


-37.4 



-74.2 

-74.2 

-74.2 




0.0107 

0.0112 

0.0120 


-74.3 


— 


0.0127 


-74.5 

-74.3 




0.0136 

0.0147 


-73.9 


— 


0.0158 


-73.0 

-71.1 

-68.5 




0.0172 

0.0188 

0.0209 


-65.3 


— 


0.0236 


-62.1 


— 


0.0267 


-59.3 


— 


0.0299 


-57.0 


-76.0 


0.0331 


-54.7 


-71.8 


0.0355 


-54.0 


-70.5 


0.0324 


-54.1 


-71.3 


0.0282 


-54.7 


-70.2 


0.0265 


-54.7 


-64.8 


0.0310 


-55.3 


-63.0 


0.0387 


-57.4 


-67.2 


0.0501 


-58.7 


- 72.3 


0.0631 


-58.0 


-75.5 


0.0794 


-57.5 


- 79.0 


0.0832 


-58.5 


— 


0.0741 


059.9 


— 


0.0596 


-63.5 


— 


0.0480 



_ 3 ; 59 

-2.93 

-2.59 

-1.60 

0.0010 -1.27 

0.0013 -0.94 

0.0019 -0.61 

0.0028 -0.28 

0.0039 0.05 

0.0050 0.38 

0.0065 0.0009 0.71 

0.0071 0.0011 1.04 

0.0070 0.0010 1.37 

0.0066 0.0011 1.70 

0.0065 0.0020 2.03 

0.0061 0.0025 2.36 

0.0048 0.0010 2.69 

0.0041 0.0009 3.02 

0.0045 3.35 

0.0047 3.68 

0.0042 4.02 

0.0036 4.35 

0.0024 4.68 
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1126.70 
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030 


6 
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349.34 
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100 


9 
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1135.27 


1137.58 


247.14 
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TABLE VIII (Fig. 15) 



mode 


tine 


f 

u 


f L 


f 

r 


0 


E 


010 


0 


561.20 


558.27 


559.73 


190.78 
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020 
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1120.66 


1122.96 


244.76 
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15 
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25 

26 

27 

28 

29 

30 

31 

32 

33 
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35 
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37 
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553.7 


-52.7 


-77.8 




0.0082 


0.0005 




-4.17 






554.2 


-52.4 


- 78.2 


0.0086 




-3.83 








554.7 


-51.8 


-78.2 


0.0091 


-3.49 








555.2 


-51.4 


-78.4 


0.0095 


-3.15 








555.7 


-51.0 


-79.0 


0.0100 


-2.81 








556.2 


-50.5 


-79.5 


0.0107 


-2.47 








556.7 


-49.7 


- 79.5 


0.0116 


-2.13 








557.2 


-49.1 


-78.9 


0.0124 


-1.79 








558.2 


-47.8 


-75.7 


0.0145 


-1.12 








558.7 


-47.0 


-72.7 


0.0158 


-0.78 




0.0011 




559.2 


-46.0 


-70.0 


0.0178 


-0.44 






559.7 


-45.0 


-67.3 


0.0200 


0.0015 


-0.10 






560.2 


-43.6 


-65.3 


0.0234 


0.0019 


0.24 






560.7 


-42.6 


-64.0 


0.0264 


0.0022 


0.58 






561.2 


-41.0 


-61.7 


0.0316 


0.0029 


0.92 






561.7 


-40.5 


-69.4 


0.0335 


0.0012 


1.25 






562.2 


-41.0 


-68.5 
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0.0013 


1.59 






562.7 


-43.2 


-58.3 


0.0247 


0.0043 
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563.2 


-46.0 
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0.0178 


0.0041 


2.27 
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-49.0 
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0.0013 
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0.0094 
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-52.4 
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-61.9 
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0.0085 


0.0029 
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0.0105 


0.0105 


0.0035 


4.64 


567.2 
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0.0056 
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